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Low-lying hypernuclei in the relativistic quark-gluon model 
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Low-lying hypernuclei \H, %H, \He, %He are described by the relativistic nine-quark equations 
in the framework of the dispersion relation technique. The approximate solutions of these equations 
using the method based on the extraction of leading singularities of the amplitudes are obtained. 
The relativistic nine-quark amplitudes of hypernuclei, including the quarks of three flavors (u, d, 
s) are calculated. The poles of these amplitudes determine the masses of hypernuclei. The mass of 
state \H with the isospin 1 = and the spin-parity J p = i + is equal to M = 2952 MeV. 
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O ■ I. INTRODUCTION. 

Hypernuclei spectroscopy is enjoying an experimental renaissance with ongoing and planned program at DA$NE, 
FAIR, Jefferson Lab, J-PARC and Mainz providing motivation for enhanced theoretical efforts (for a recent review, 
seeRef. Q). 

LQCD calculations are performed using an isotropic quark action at the SU (3) flavor symmetric point corre- 
sponding to the physical strange quark mass, with — vhk — fn v ~ 800 MeV. The lattice spacing b = 0.145 fm 
has been used in the calculations, dictated by the avaible computational resources, and therefore an extrapolation to 
the continuum has not been performed. Further, extrapolation of the physical pion mass have not been attempted, 
because the quark mass dependence of the energy levels in the light nuclei are not known. Future calculations at 
smaller lattice spacing and at lighter quark masses will facilitate such extrapolations and lead to first predictions for 
the spectrum of light nuclei, with completely quantified uncertainties, that can be compared with experiment. 

In the recent paper Q the relativistic six-quark equations are found in the framework of coupled-channel formalism. 
The dynamical mixing between the subamplitudes of hexaquark are considered. The six-quark amplitudes of dibaryons 
are calculated. The poles of these amplitudes determine the masses of dibaryons. We calculated the contribution of 
six-quark subamplitudes to the hexaquark amplitudes. 

In the previous paper Q the 3 He as the system of interacting quarks and gluons is considered. The relativistic 
nine-quark equations are found in the framework of the dispersion relation technique. The dynamical mixing between 
the subamplitudes of 3 He is taken into account. The relativistic nine-quark amplitudes of 3 He, including the u, d 
\ quarks are calculated. The approximate solutions of these equations using the method based on the extraction of 
leading singularities of the amplitudes were obtained. The poles of the nonaquark amplitudes determined the mass 
of 3 He. The mass of 3 He state is equal to M — 2809 MeV. The model use only two parameters: the cutoff A = 11 
[9l [To| and gluon coupling constant g = 0.1536, which is smaller than the gluon coupling constant of baryon physics. 

In the present paper the hypernuclei \H , %H , \He, ^He, nnA and nnS are considered in the framework of the 
dispersion relation technique. The approximate solutions of the nine-quark equations using the method based on 
the extraction of leading singularities of the amplitude are obtained. The relativistic nona-amplitudes of low-lying 
hypernuclei, including the three flavors (u, d, s) are calculated. The poles of these amplitudes determine the masses of 
the hypernuclei. In Sec. II the relativistic nine-quark amplitudes of the hypernuclei are constructed. The dynamical 
mixing between the subamplitudes of hypernuclei are considered. 

Sec. Ill is devoted to the calculation results for the masses of the lowest hypernuclei (Table HJ). In conclusion, the 
status of the considered model is discussed. 



'Electronic address: gerasyuta@SG6488.spb.edu 
^Electronic address: matskev@pobox.spbu.ru 



2 



TABLE I: Masses of lowest hypernuclei. 
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II. NINE-QUARK AMPLITUDES OF HYPERNUCLEI. 



We derive the relativistic nine-quark equations in the framework of the dispersion relation technique. We use only 
planar diagrams; the other diagrams due to the rules of 1/N C expansion [llTfl3l ] are neglected. 

The current generates a nine-quark system. The correct equations for the amplitude are obtained by taking into 
account all possible subamplitudes. Then one should represent a nine-particle amplitude as a sum of 36 subamplitudes: 



i<3 



This defines the division of the diagrams into groups according to the certain pair interaction of particles. The total 
amplitude can be represented graphically as a sum of diagrams. We need to consider only one group of diagrams and 
the amplitude corresponding to them, for example A 12 . We consider the derivation of the relativistic generalization 
of the Faddeev-Yakubovsky approach. 

In our case the two states of hypertritons are obtained: with the isospin 1 = (uud udd uds) and with the 
isospin 1=1 (uududduus). The similar method allows us to construct the hyperhelium states: with the isospin 
1 = %fle (uuduuddds), with the isospin 1=1 \He (uuduududs) and with the isospin 1 = 2 fjiJe (uud uud uus) . 

We take into account the pairwise interaction of all nine quarks. The set of diagrams associated with the amplitude 
A±2 can further be broken down into some groups corresponding to subamplitudes similar to [8|, Appendix A. 

In order to represent the subamplitudes A+ in the form of a dispersion relation, it is necessary to define the amplitude 
of qq interactions. We use the results of our relativistic quark model [H |: 

I \ G n (Sik) , . 

a n (s ik ) = - -—, r, (2) 

1 - Bn(Sik) 



B n ( Slk) = [ d Jkpj^m, (3 ) 

J 7T S ik - S lk 

(mi+m k ) 2 

Here G n (sik) are the diquark vertex functions ( Table ILT| . The vertex functions are determined by the contribution of 
the crossing channels. The vertex functions satisfy the Fierz relations. These vertex functions are generated from gy . 
B n (sik) and p n (sik) are the Chew-Mandelstam functions with cutoff A and the phase spaces: 

Pn ( Slk ,J P ) = U(n,J p ) ( ^ ^ +P(n,J p ) + S(n,J p ) {m *- mk)2 
\ (m l + m k )' L Sik 



x \J (sik - (rrii + m k ) 2 )(s ik - (m, - m k ) 2 ) 



The coefficients a(n, J p ), /3(n, J p ) and 5(n, J p ) are given in Table [III 



3 



TABLE II: The vertex functions and coefficients of Chew-Mandelstam functions. 
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Here n = 1 coresponds to gg-pairs with J p = + , n = 2 corresponds to the pairs with J p = 1+. is the 
two-particle subenergy squared. 

In the case in question the interacting quarks do not produce a bound states, thefore the integration is carried out 
from the threshold (m, + m^) 2 to the cutoff A. 

Let us extract singularities in the coupled equations and obtain the reduced amplitudes a*. 

The isospin 1 = \H. 

We need to consider the system of 29 equations: 



ai : l uu , l M , 0" d , 0" s , ds (5) 

£^ . -^uu-^uu -^uu-^dd -^uiiQud -^uuqus ^uu^ds 
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The isospin I = 1 |if (24 equations) 



ai : 1™, l dd , ud , US , ds (8) 
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Here the ai are determined by the diquarks, the c%2 includes the two diquarks and the five quarks, as corresponds 
to the ppA (pnA) states. 

The equations for the \He, ^He are considered. 

The isospin 1 = 0. 

We need to consider 23 equations. 
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The isospin 7 = 1 (25 equations): 



a x : l uu , r d , ud , US , ds (14) 
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The isospin 1 = 2 (20 equations): 



ai : l uu , l dd , ud , US , ds (17) 
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The coefficients of the coupled equations are determined by the permutation of quarks (l5l . Il6j . For the simplicity 
we consider the graphical equation of the reduced amplitude a\ 1 1 (Fig- 1)- 

In Fig. 1 the first coefficient is equal to 12, that the number 12 = 2 (permutation of particles 1 and 2) x2 
(permutation of particles 3 and 4) x3 (permutation of pairs (56) and (12), (56) and (34)); the second coefficient equal 
to 12, that the number 12 = 2 (permutation of particles 1 and 2) x3 (permutation of pairs (12) and (34), (12) and 
(56)) x2 (we can replace 7-th d-quark with 8-th <i-quark); the third coefficient equal to 6, that the number 6 = 2 
(permutation of particles 1 and 2) x3 (permutation of pairs (12) and (34), (12) and (56)); the 4-th coefficient equal to 
6: 6 = 2 (permutation of d-quarks 7 and 8) x3 (permutation of pairs (12) and (34), (12) and (56)); the 5-th coefficient 
equal to 12: 12 = 2 (permutation of ci-quarks 7 and 9) x2 (permutation of d-quarks 7 and 8) x3 (permutation of 
pairs (12) and (34), (12) and (56)); the 6-th coefficient equal to 24: 24 = 2 (permutation of particles 1 and 2) x2 
(permutation of particles 3 and 4) x3 (permutation of pairs (56) and (12), (56) and (34)) x2 (permutation of d-quarks 
7 and 8); the 7-th coefficient equal to 48: 48 = 2 (permutation of particles 1 and 2) x2 (permutation of particles 3 and 
4) x3 (permutation of pairs (56) and (12), (56) and (34)) x2 (permutation of (i-quarks 7 and 9) x2 (permutation of 
d-quark 7 and s-quark 8); the 8-th coefficient equal to 48: 48 = 2 (permutation of particles 1 and 2) x2 (permutation 
of particles 3 and 4) x2 (permutation of pairs (12) and (34)) x3 (permutation of pairs (56) and (12), (56) and (34)) 
x2 (permutation of d-quarks 7 and 8); the 9-th coefficient equal to 24: 24 = 2 (permutation of particles 1 and 2) x2 
(permutation of particles 3 and 4) x2 (permutation of pairs (12) and (34)) x3 (permutation of pairs (56) and (12), 
(56) and (34)); the 10-th coefficient equal to 24 (this case is similar to 6-th); the 11-th coefficient equal to 48 (this 
case is similar to 7-th); the 12-th coefficient equal to 48: 48 = 2 (permutation of particles 1 and 2) x6 (permutation 
of pairs (12), (34) and (56)) x2 (permutation of d-quarks 7 and 8) x2 (permutation of d-quark 8 and s-quark 9); the 
13-th coefficient equal to 24: 24 = 2 (permutation of particles 1 and 2) x6 (permutation of pairs (12), (34) and (56)) 
x2 (permutation of d-quarks 8 and 9). 

The similar approach allows us to take into account the coefficients in all equations. 

The functions I\ - Iio are taken from the paper 0]. The other functions are determined by the following formulae 
Q (see also Fig. 1): 

7ii(12,34,15,36,47) = 7i(12,15) x 7 2 (34, 36, 47) , (20) 



7 12 (12,34,56,17) = 7 1 (12,17), (21) 
713(12, 34, 56, 17, 28) = 7 2 (12, 17, 28) , (22) 
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7^(12, 34, 56, 17, 38) = 7 X (12, 17) x 7^34,38) , (23) 

7i 5 (12, 34, 56, 23, 47) = J 7 (12, 34, 23, 47) , (24) 

7i 6 (12, 34, 56, 17, 23, 48) = 7 8 (12, 34, 17, 23, 48) , (25) 

ii 7 (12,34,56,17,45) = A (12, 17) x 7 3 (34, 56, 45) , (26) 

lis (12, 34, 56, 17, 38, 49) = h (12, 17) x 7 2 (34, 38, 49) , (27) 

7i 9 (12,34,56,17,38,59) = /i(12,17) x 7i(34,38) x 7i(56,59), (28) 

7 20 (12, 34, 56, 17, 45, 68) = h (12, 17) x 7 7 (34, 56, 45, 68) , (29) 

hi (12, 34, 56, 17, 28, 45) = 7 2 (12, 17, 28) x 7 3 (34, 56, 45) . (30) 
The main contributions are determined by the functions I\ and 7 2 : 
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tum - B MVBk(s?) f ds' 12 Gi(s^) Pi (s[ 2 ) 1 fd Zl (2) f 
2( J ' ~ B t (s?) J n s' 12 -stf 2nJ 2 J 
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y/1 - z((2) - z 2 (2) - ^ 2 (2) + 2z!(2)z 2 (2)z 3 (2) 



1- B,( S ' 13 ) 1-B k (s' 24 )> 



(32) 



where i, j, k correspond to the diquarks with the spin-parity J p = + , 1 + . 

The one new contribution 7 22 ~ 10 -8 . The contributions of Iiq, In, 7ig, I 2 o, I 2 \ also are small ~ 10~ 7 . We do not 
take into account these functions in the coupled equations. 

III. CALCULATION RESULTS. 

The poles of the reduced amplitudes correspond to the bound state and determine the masses of low-lying hyper- 
nuclei. The quark masses m = 410 MeV and m s = 557 MeV coincide with the quark masses of the ordinary baryon 
in our model [H, [l(| ■ This model does not have a new parameter as compared with the Ref. Q . The gluon coupling 
constant g = 0.1536 is determined by fixing the mass of 3 i7e at M = 2809 MeV. The cutoff parameter is similar to 
ordinary value A = 11. The calculated mass values of the hypernuclei are shown in the Table [I] 

We predict the degeneracy of hypertriton and hyper helium with the spin-parity J p = i + , | + . 

The hypertriton with the 7 = \H and the hyperhelium with the 7=1 \He also can be considered with the 
degeneracy. 
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IV. CONCLUSION. 



It is now clear, but hardly a surprise, that the spectrum of nuclei and hypernuclei changes dramatically from light 



quark masses. While one had already shown this from the recent work [17| on the _ff-dibaryon, and nucleon-nucleon 
scattering lengths, this has now been demonstrated to the true for even larger systems. 

It will be interesting to learn how the various thresholds for binding evolve with the light quark masses providing 
accurate binding energies for any given light quark masses will require the inclusion of electromagnetic effect. 

Lattice QCD has evolved to the point where first principles calculation of light nuclei are now possible, as demon- 
strated by the calculation at unphysically heavy quark masses presented in the paper [lS]. The experimental program 
in hypernuclear physics, and the difficulties encountered in accurately determined rates for low energy nuclear reac- 
tions, warrant continued effort in, and development of, the application of LQCD to nuclear physics. 

Cleary, calculation at smaller lattice spacings at SU (3) symmetric point are required in order to remove the 
systematic uncertainties in the nuclear binding energies at these quark masses. In order to impact directly the 
experimental program in nuclear and hypernuclear physics, analogous calculations must be performed at lighter 
quark masses, ideally at or close to their physical values. 

The antihypertriton have the mass that is similar to the hypertriton in our calculations. 

The experimental data of hypertriton mass is about 3 GeV [l9| . 
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